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Abstract 

We review theoretical calculations for top-quark production that include complete 
next-to-leading-order QCD corrections as well as higher-order soft-gluon corrections from 
threshold resummation. We discuss in detail the differences between various approaches 
that have appeared in the literature and review results for top-quark total cross sections 
and differential distributions at the Tevatron and the LHC. 



1 Introduction 



Top-quark physics is a centerpiece of the research programs at the Tevatron and the LHC. 
Due to its unique position as the heaviest particle known to date, the top quark not only plays 
a role in many models of physics beyond the Standard Model (SM), but also in the precision 
electroweak fits constraining the mass of the Higgs boson. 

The dominant mechanism for creating top quarks at hadron colliders is the production 
of a top-antitop pair through QCD interactions. Tens of thousands of top-quark pairs have 
already been produced and studied at the Tevatron since the discovery of the top quark there 
in 1995 [U and millions will be produced at the LHC. As the measurements become more 
precise, it will become increasingly important to have reliable QCD predictions of the total 
and differential pair-production cross sections. This is true not only if signals of new physics 
manifest as slight excesses in differential cross sections, but also if top-quark properties are 
completely determined through the SM. In that case top-quark pair production will be a 
benchmark process, used in tasks from subtracting backgrounds to constraining gluon PDFs 
in regions of x relevant for Higgs production. 

Another process, with a smaller cross section than for top-antitop, is single top-quark pro- 
duction, which can proceed via three distinct partonic channels. Single top-quark production 
is important in probing electroweak theory and studying the electroweak properties of the top 
quark as well as for the discovery of new physics, since the top-quark mass is of the same order 
of magnitude as the electroweak symmetry breaking scale. 

The main purpose of this review is to summarize certain aspects of QCD calculations for 
the total and differential pair-production cross sections. We will focus on inclusive observables 
VPip) ~^ tiX, summed over spins, although some results for the single top-quark total cross 
section in the t-channel, s-channel, and in tW production are also discussed. In particular, we 
will perform detailed discussions of the total inclusive production cross section, the top-pair 
invariant mass distribution, the transverse momentum (px) distribution of the top (or antitop) 
quark, and the rapidity distribution, as well as the forward-backward (FB) asymmetry at the 
Tevatron. While this leaves out a number of interesting topics related to other exclusive 
observables, these are covered in excellent reviews of top-quark physics in the literature, e.g. 
[3]. In this review we present in detail the most up-to-date theoretical predictions for the 
observables mentioned above. 

The starting point for a study of inclusive observables in top-quark pair production is the 
next-to-leading-order (NLO) calculation of differential and total cross sections performed more 
the two decades ago [U El [7] . However, to keep up with experimental precision requires to 
go beyond them. Many theorists have responded to this challenge by tackling parts of the 
diagrammatic calculations needed to reach next-to-next-to-leading-order (NNLO) accuracy 
and a large number of papers have appeared with pieces of the NNLO calculation. Due to 
very recent progress it now looks feasible that the NNLO cross section may be available in the 
near future, a step forward that will mark a major accomplishment. 

In the absence of the full NNLO results, an important tool for including higher-order cor- 
rections is soft gluon resummation at next-to-leading-logarithm (NLL) accuracy [HI E] and 
beyond. The resummation techniques apply to both differential and total cross sections and 
fixed-order expansions were derived from the NLL resummed cross section with additional 
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subleading logarithms in [IU1EEI]. Recently, next-to-next-to-leading-logarithm (NNLL) predic- 
tions have appeared for the total cross section [T2J [131 El US] and for double differential cross 
sections with respect to the top quark pr and rapidity [HI [151 HB] , or the pair invariant mass 
and rapidity [13]. A primary goal of this review is to clarify the similarities and differences 
between the different approaches. 

The remainder of this document is organized as follows. In Section [2] we review the basics 
of fixed-order calculations to NLO and summarize progress up to NNLO. In Section |3l we 
discuss the underlying assumptions and techniques in soft gluon resummation, providing a 
theory comparison of different approaches. We present numerical results for the total inclusive 
cross section within various approaches in Section [H for a selected group of differential cross 
sections in Section \5[ and for the FB asymmetry in Section El Single top production is briefly 
discussed in Section [71 and some closing remarks are made in Section [HJ 



2 Kinematics and fixed-order results 

In this section we discuss the basics of fixed-order calculations for inclusive top-quark pair 
production at hadron colliders. We thus consider the scattering process 

Nx(Pi) + N 2 (P 2 ) -> i(p 3 ) + t(p 4 ) +X, (1) 

where N± and N 2 indicate the incoming protons (LHC) or proton and anti-proton (Tevatron), 
and X represents an inclusive hadronic final state. The calculation of the differential cross 
section relies on factorization in QCD. This is the statement that a generic differential hadronic 
cross section is given in terms of partonic cross sections associated with scattering of gluons and 
quarks, after a convolution integral with parton distribution functions (PDFs) that describe 
the parton content of the hadrons. We write this schematically as 

dcr = ^2 J dxidx 2 ^ i / Nl {x 1 ,n F ) (f} j / N2 (x 2 ,fi F ) dfr i j(x 1 ,x 2 ,[J,F,lJ,R,a s (iJ, R )), (2) 

where the Xi are longitudinal momentum fractions of the incoming partons i, <fii/N are PDFs, 
and the d&ij are partonic cross sections. The collinear singularities are factorized in a process- 
independent manner and absorbed into the PDFs, which depend on the factorization scale 
fi F . The partonic cross sections can be expanded in a fixed-order series in the strong coupling 
constant a s (^R) as 




where the first term in square brackets is referred to as leading-order (LO), the second term 
next-to-leading-order (NLO), the third term next-to- next-to-leading-order (NNLO), and so on. 
The physical cross section is independent of the factorization scale /ip and the renormalization 
scale fiR. However, the truncation of the infinite perturbative series at finite order typically 
results in a non-negligible numerical dependence. 
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Figure 1: The four Feynman diagrams contributing to top-quark pair production at LO. 



At LO, the partonic cross sections receive contributions from the quark- ant iquark annihi- 
lation and gluon fusion subprocesses 

q(pi) + q(p 2 ) -> t(p 3 ) + f(p 4 ) , 

g(pi) + gfa) -> t( P3 ) + i( P4 ) . (4) 

The Feynman diagrams needed for the computation of the differential cross section are shown 
in Figure [TJ 

The results for the LO differential partonic cross sections read 



~Jq°g txC f ( i\ + u\ 2m, 2 . 
dtidui N c \ s 2 s 



S dkdu X (iV2 - 1)2 ^ F t l u 1 



S 2 S tiUi 



Sis + h + ut), (5) 



where C F = (iV c 2 - 1)/2N C and C A = N c , with N c = 3 the number of colors in QCD. The 
partonic Mandelstam variables are defined as 

s = Oi+p 2 ) 2 , k = (p! - p 3 ) 2 - m\ , u x = (P2 - Psf ~ m 2 t , (6) 

where m t is the top quark mass. The partonic momenta are related to the hadronic ones 
through pi = X\P\ and pi = X2P2, so the connection between the hadronic invariants 

s = (P l + P 2 ) 2 1 tl = (P l -p 3 ) 2 -m 2 t , Ul = (P 2 -p 3 ) 2 -m 2 tJ (7) 

and the partonic ones is 

s = X1X2S , t\ = xiti , u\ = X2U1 . (8) 

At Born level, any differential cross section can be derived from (j5J) by an appropriate 
change of variables. For instance, single-particle-inclusive (1PI) cross sections depending on 
the rapidity y and px of the top quark in the laboratory frame are calculated using 

t\ = — y/s m_i_ e~ y , u\ = — y/sm±e y , (9) 

where m± = \Jp\ + m 2 , while pair-invariant-mass (PIM) observables depending on the pair- 
invariant mass M t | = (p 3 + P4) 2 and scattering angle of the top-quark are obtained using 

M 2 - M 2 - 
s = M 2 t , tl = —£-(l-p t ccxie), U! = — -*(1 + Acos0). (10) 
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where (3 t = a/1 - 4rrij/M%. 

The total cross cross section itself is conventionally written in the form 



a(s,m t ) 



in 



t ij - </4m 



4mj 



where the parton luminosities are defined as 

^ dx 



f dx 

®ij{y,VF)= / — 4>i/ Nl (x,fi F )(f) j/N2 (y/x,fi F ) . 

J V x 



(12) 



The perturbative scaling functions are easily obtained by integrating the differential par- 
tonic cross sections over the appropriate phase space. We define their fixed-order expansion 

as 



f« = 4° 3 + l-nj;/ + (1-n J 2 + . . . . 
These scaling functions are dimensionless quantities normally written in terms of 



Am 2 



P 



(13) 



(14) 



where at Born level /3 is just the magnitude of the three-velocity of the top or antitop quark. 
For reference, the LO expressions in the two channels read 



(0) _ ixN c C F (3p / 1 



(15) 



C f 



P 2 \ „ P 2 ' 
P ~2J + A 4 



In 



^|)-^(l + p)-§(4 + 5p) 



(16) 



Predictions made using the LO partonic cross sections depend very strongly on the choice 
of the factorization and renormalization scales and are not appropriate for any sort of detailed 
phenomenology. However, they are sufficient for understanding most of the qualitative features 
of the observables discussed later on. For instance, the LO results predict a total cross section 
of roughly 7 pb at the Tevatron, with a dominant contribution of around 90% from the qq 
channel. At the LHC with y/s = 7 TeV the cross section is around twenty times as large, 
and the dominant contribution, of around 75%, is from the gg channel. Obviously, the size of 
the cross section and relative strength of the different channels reflects the differences in the 
parton luminosities at the two colliders, since the partonic cross sections are the same. 

To make quantitative comparisons between theory and experiment requires the partonic 
cross sections at NLO accuracy. The total and differential inclusive production cross sections 
were calculated at NLO more than 20 years ago [U El El [7]. All of the NLO results for 
top-quark pair production are implemented in parton Monte Carlo programs such as MCFM 
[TT] or the NLO version of MadGraph/MadEvent based on [TSl [T9"] . or also including parton 
showers in programs such as MC@NLO [2U]. These programs allow for the calculation of 
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arbitrary differential distributions and are thus an important tool for experimentalists and 
theorists involved in phenomenological applications. 

Numerical results using the partonic cross sections at NLO will be reviewed in the follow- 
ing sections. It is conventional to estimate the theoretical uncertainties associated with the 
truncation of the perturbative series to NLO by picking default values of the factorization and 
renormalization scales, varying them up and down by factors of two, and taking the resulting 
spread of values for the cross section as an indication of the size of uncalculated, higher order 
terms. Although significantly more stable than LO predictions under such scale variations, 
NLO calculations still suffer from scale uncertainties of roughly 10-20%. Higher-order calcu- 
lations using soft-gluon resummation reduce the scale uncertainty and will be covered in the 
next section. An even more ambitious approach is to calculate the full NNLO corrections, and 
we now review recent progress in this area. 

The calculation of the NNLO cross sections involves two sets of corrections: i) virtual 
corrections, which can be split into genuine two-loop or one-loop squared contributions, and 
ii) real radiation, which involves either one- loop diagrams with the emission of one extra parton 
in the final state, or tree-level diagrams with two extra partons in the final state. Needless to 
say, the NNLO calculation is very challenging, and so far only pieces of it are available. 

For the virtual corrections, the current situation can be summarized as follows. Analytic 
results for the two- loop contributions in the high-energy limit are known from [2H 122] , along 
with contributions to the qq channel from fermionic [23] and planar corrections [23] , and to the 
gg channel from planar corrections [251 . The exact results for the two- loop contributions in 
the qq channel are known numerically from [2S] • In addition, the two-loop infrared singularity 
structure in the qq and gg channels was determined in [23 EH]. Finally, one-loop squared 
terms were calculated in [211 EOj [31] . 

For the real corrections, the pieces arising from one-loop virtual-real graphs are known 
from calculations of tt production in association with jets [321 13"3"| 13"^ I3"B] . The more involved 
task of developing a new subtraction scheme and calculating the contributions from double 
real radiation has been dealt with in [361 EZl EH EH1 SO]. In light of these advances, the 
completion of the NNLO calculation now looks feasible, although a significant amount of work 
is still required to assemble all the elements. In the meantime, approximations of the NNLO 
corrections are available using techniques from soft gluon resummation. 

We end this section with a couple of comments on progress in NLO calculations which 
is slightly outside the main stream of presentation but nonetheless very important. So far, 
we have considered only inclusive top-pair production with a final state ttX. In reality, 
the top quark decays almost as soon as it is produced, and experiment observes the decay 
products rather than the top quarks themselves. If the theory calculations stop at the level 
of on-shell top-quark production, experimentalists must fill in the gaps by correcting their 
measurements of the decay products (with a possibly complicated set of acceptances) back 
to that level. Obviously, this should be done as accurately as possible. Recent calculations 
of top-quark production at NLO, including the top-quark decay in the dominant t — > bW 
channel (along with leptonic decay of the W), allow for a much more direct comparison 
between theory and experiment. The most advanced computations [HI [32] take into account 
the full correlations between production and decay, without the further assumption of on-shell 
top-quark production used in the narrow-width approximation. A main result of those works 
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is that the narrow-width approximation is valid to the percent level for sufficiently inclusive 
observables (including those studied in this review). This means that in most cases one may 
use equally well the previous calculations of [US HI] , which include a complete description of 
production and decay to NLO within the narrow-width approximation. 

3 Soft gluon resummation 

So far we have considered only the fixed-order perturbative expansion for the partonic cross 
sections. However, there are certain cases where such an expansion is not the optimal calcu- 
lational procedure, due the presence of large logarithms related to soft gluon emission. The 
reorganization of the perturbative expansion in these cases is carried out through techniques 
referred to as soft gluon resummation. 

Soft gluon resummation was first applied to top-quark production at leading-logarithm 
(LL) accuracy twenty years ago in [15]. At that level the resummation is universal and only 
depends on the identity of the incoming partons (quarks and gluons) and does not depend 
on the details of the hard scattering. Other approaches at LL followed in [16J H7] and [H] . 
However, beyond LL the color structure of the process directly enters the resummation. Re- 
summation for top quark production at NLL was first presented in [SI IS] where the one-loop 
soft anomalous dimension matrices were calculated. NNLO expansions of the differential re- 
summed cross section, incorporating additional subleading logarithms, were presented in 1PI 
and PIM kinematics in [10] . A different approach at NLL for the total cross section only, using 
production threshold, appeared in [19], and was used in [50] and with additional subleading 
logarithms (and Coulomb terms) in |51j . 

Further work on two-loop calculations for soft and collinear singularities in massless gauge 
theories [521 lEU G31 E5J M\ |57] and with massive quarks [5S1 ESI EDI EU E21 [23 I2S1 H31 E31 EH 
HU [15] have eventually allowed NNLL resummation for ti production [121 021 HH EEl [16] . 

A typical example where resummation is needed is the top-pair invariant mass distribution 
at large values of the invariant mass M t i . To understand why this is the case, consider the 
factorized expression for the hadronic cross section, which reads 

da f 1 dz daij(z,M th ii F ^ R ,a s (ii R )) 

I yW^) dWu • (17) 

In the limit of very large invariant mass, the variable r = M^/s —> 1, which implies that 
the partonic variable z = M^/s — > 1. In that limit, the partonic center-of-mass energy is 
just above the threshold required to produce the top pair with a given invariant mass, which 
forces any additional radiated partons to be soft. Such soft radiation is associated with IR 
divergences in the QCD scattering amplitudes, and after collinear subtractions and cancel- 
lations with singularities from the virtual corrections, the cross section contains logarithmic 
plus-distribution corrections. For the n-th order correction to the partonic cross sections in 
the gg and qq channels, these are of the form 

; m = 0, ••• ,2n-l, (18) 

+ 



ln m (l 
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where the plus distributions are defined as 



dz 



ln m fl 



g(z)= I 1 dz ln "j ( ^ ~ * } W) - 9(1)] + ^ ln m+1 (l - r) (19) 



for an arbitrary function g(z). The qg channel also receives up to double logarithmic correc- 
tions but it is suppressed by a relative factor of (1 — z) compared to the gg and qq channels. 

The large logarithms appearing in the limit z — > 1 make the perturbative series poorly 
behaved, and to make reliable predictions requires that they be resummed. Because the 
logarithms are related to the fact that at the partonic threshold real gluon emission is soft, such 
a resummation is interchangeably referred to as "soft gluon" or "threshold" resummation. The 
technical machinery required to perform soft gluon resummation will be reviewed in following 
two subsections. For now, we just note that the plus-distribution corrections stem from terms 
depending on the dimensionless ratio 2E S / fi, with 2E S = s(l — z) the energy of the soft-gluon 
radiation in the partonic center-of-mass frame, so the renormalization-group is the basic tool 
for resummation. 

While interesting from the technical point of view, the situation described above is of little 
phenomenological importance. The reason is that the differential cross section is essentially 
zero in the limit r — > 1, since in that case the combined energy of the top-quark pair approaches 
the collider energy and the probability that two initial state partons carry such a large fraction 
of the energy is tiny. Both at the Tevatron and at the LHC, the invariant mass distribution is 
largest in regions of phase space where r is closer to zero than to one, and the partonic cross 
section in the convolution integral (117)1 is evaluated at values of z far from unity. In that case, 
the leading terms in the soft limit give the largest contributions to the hadronic cross section 
if the parton luminosities as a function of r/z are by far largest at z — > 1. If that is true, 
then the expansion of the partonic cross section in the z — > 1 limit under the integral ( [T7J is 
parametrically justified. To a good degree this is an actual property of the parton luminosities, 
and detailed studies show that, at least at NLO, the leading terms in the z — > 1 limit provide 
an excellent approximation to the full result. Assuming that also beyond NLO the leading 
terms in the soft limit account for the bulk of the corrections, then predictions using soft gluon 
resummation are an improvement on the fixed-order expansion even in regions of phase space 
where r is not close to unity. 

The line of reasoning used for the invariant mass distribution can be applied to soft gluon 
resummation for other differential partonic cross sections, or else to the total partonic cross 
section directly. In fact, the vast literature on soft gluon resummation in top-quark pair 
production can be broken down into the three main cases shown in Table [TJ In each case, 
one considers the soft limit of the (differential) partonic cross section shown in the table, 
resums corrections which become large in that limit, and relies on the fall-off of the parton 
distributions away from the soft limit to dynamically enhance the partonic threshold region. 
The cases of PIM and 1PI kinematics were first considered in [H [9] and [65], respectively, 
and work at the level of double differential cross sections. The logarithmic corrections in PIM 
kinematics are of the form ( 1T9"|) . while in 1PI kinematics they are of form 



a 



ln m (g 4 /m f 2 ) 

s 4 



m = 0," ■ ,2n- 1, (20) 
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Name 


Observable 


Soft limit 


pair-invariant-mass (PIM) 
single-particle- inclusive (1PI) 
production threshold 


da/dMudO 
da/dpxdy 
a 


(!-«) = 1 - Afg/a -> 
s 4 = s + ti +ui —> 
P = y/l- 4mf/s ->> 



Table 1: The three cases in which soft gluon resummation has been applied. The first column 
indicates the name often used in the literature, the second the observable to which it applies, 
and the third the partonic variable associated with large logarithmic corrections in the soft 
limit. 

where s 4 = s + 1\ + iii and the plus distributions are defined as 

(21) 

It is very important to emphasize that resummed results in 1PI and PIM kinematics are ob- 
tained after integrating over a specific portion of the fully differential phase space with respect 
to additional soft gluon radiation. This means that the results in 1PI or PIM kinematics 
apply only to the differential cross sections in the table, in other words they are not related 
by a simple change of variables. However, they can both be used to predict the total cross 
section by integrating over the distributions. Alternatively, one can perform resummation for 
the total partonic cross section directly working in the production threshold limit, which is 
the third entry in the table. In that case the large corrections in the soft limit involve simple 
logarithms of /3 instead of singular plus distributions. 

As far as soft gluon resummation is concerned, results in the production threshold limit 
are actually a special case of PIM and 1PI kinematics and thus do not contain independent 
information. For instance, the limit (1 — z) — > plus the additional phase-space restriction 
M t i — > 2mt, or s 4 — > plus the restriction iii + t\ = —4mf, both imply the limit /3 — > 0. 
However, in the limit /3 — > there are additional Coulomb singularities of the form ln m f3//3 n , 
not all of which are determined by soft gluon resummation alone. In comparing results for 
the total cross section obtained in the different limits, an important consideration is whether 
the subleading terms in /3 contained in the 1PI and PIM results, or the Coulomb singularities 
are more important. This is a numerical question which we will come back to later on. In 
the meantime, we pause to explain in more detail the technical formalism needed to perform 
resummation within the different types of kinematics. 

3.1 Mellin-space resummation 

The resummation of the threshold logarithms has traditionally been performed in Mellin 
moment space. By taking moments, divergent distributions in 1 — z (PIM kinematics) or s 4 
(1PI kinematics) produce powers of In N, with iV the moment variable: 

t dz z N ~ 1 
Jo 



m m n 



m 



— ln m+1 iV + O (In™" 1 N) 



(22) 



8 



We define moments of the partonic cross section by 0"(iV) = j dz z N ~ 1 a(z) (PIM) or by 
er(iV) = j{dsn/s) e~ Ns4 / s a(s4) (1PI). Then the logarithms of N in cr(N) exponentiate. 

Consider the partonic process /1 + /2 — > t+X where X represents the additional final-state 
particles apart from a produced top quark. We proceed with the derivation of the resummed 
cross section by writing a factorized form for the moment-space parton-parton scattering cross 
section, infrared-regularized by e, <j/ 1 / 2 ^ i x(A r , e), which factorizes as the hadronic cross section 

°hh-+tx(N, e) = 4> fl / fl (N, n F , e) 4>f 2 /f 2 (N, fi F , e) a flh ^ tx {N, fi F , fj, R ) , (23) 

where the moments of <p are given by </>(iV) = J Q dx x N ~ 1 (f)(x), and for simplicity we do not 
show the dependence of the cross section on kinematical variables. We factorize the initial- 
state collinear divergences into the parton distribution functions, </>, which are expanded to 
the same order in a s as the partonic cross section, and thus obtain the perturbative expansion 
for the infrared-safe partonic short-distance function a. 

The partonic function a is still sensitive to soft-gluon dynamics through its N dependence. 
We then refactorize the moments of the cross section as [H] 

(T hh ^ tx {N, e) = $ fl/fl (N, up, e) vj> h/h (N, fi F , e) 
x H^ tx ( a .fa)) S f L Y^ tx (J^-., a.Qi R )j ]J J 3 (N, p F , e) + 0{l/N) , (24) 

where i/j are center-of-mass distributions that absorb the universal collinear singularities from 
the incoming partons, Hjl is an N- independent function describing the hard-scattering, Su 
is a soft gluon function associated with non-collinear soft gluons, and J are functions that 
absorb the collinear singularities from any massless partons in the final state. We note that 
the J functions do not appear in the resummation for tt production but are needed in single 
top production via the t- and s-channels. 

H and S are matrices in the space of the color structure of the hard scattering, with color 
indices / and L. The hard-scattering function involves contributions from the amplitude of 
the process and its complex conjugate, Hn = h* L hi. 

Using Eqs. ( 123]) and (124)) we can write a in terms of H, S, J, and the ratios ip/(f). 
The constraint that the product of these functions must be independent of the gauge and 
factorization scale results in the exponentiation of logarithms of N in the parton and soft 
functions EE] . 

The soft matrix Sli depends on N through the ratio M/(Nfi F ), and it requires renormal- 
ization as a composite operator. However, in the product Ha Sli the UV divergences of S 
are balanced by those of H. Sli satisfies the renormalization group equation [U [9] 

H— + P(g.)— Sli = -(^ s ) LB S B i ~ Sla(Ts)ai , (25) 



dfx dg 

where (3(g s ) is the QCD beta function and g 2 s = 47ia s . Ts is the soft anomalous dimension 
matrix, and it is calculated in the eikonal approximation by explicit renormalization of the 
soft function. In a minimal subtraction renormalization scheme with e = 4 — n, Ts is given at 
one loop by 

T { s\9s) = -^£-Res^ Q Z s (g s ,e) . (26) 
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In processes with simple color structure is a 1 x 1 matrix while in processes with complex 
color structure it is a non-trivial matrix in color exchange. For quark- antiquark scattering 
into a top-antitop pair, is a 2 x 2 matrix [3 E| ; for gluon-gluon fusion into a top-antitop 
pair it is a 3 x 3 matrix [9]. 

The exponentiation of logarithms of N in ip/4> and J together with the solution of the renor- 
malization group equation (125 p . provide us with the complete expression for the resummed 
partonic cross section in moment space [9] 



a res (N) 



cxp 



i=l,2 



cxp 



exp 



E 2 

i=l,2 j Vf 



x 



Tr | H^ tx (a s (VI)) exp ^ V 



d V r t/x 



>.(//)) 



x exp 



VS/N' 



(27) 



In 1PI kinematics Ni = N(—ti/M 2 ) for incoming partons i, and iV' = N(s/M 2 ); here M is 
any chosen hard scale relevant to the process, such as the top quark mass. In PIM kinematics 
Ni = N' = N. Also note that N = Ne lE , with ^ E the Euler constant. 

The first exponent in Eq. (1271) arises from the exponentiation of logarithms of N in the 
ratios ip/<f) (the sum % — 1, 2 is over incoming partons), and it is given in the MS scheme by 



dz- 



^ dX 
A 



Ai(a s (Xs)) + Di [a s ((l 




(28 



where A EB] and D [HE] have well-known perturbative expansions. 

The second exponent in Eq. (12 7j) arises from the exponentiation of logarithms of iV in the 
functions Jj for final-state particles (the sum j is over outgoing quarks and gluons) [671 ESj • 
It does not appear in ti production but is needed for single top quark production. 

The third exponent in Eq. (|27p controls the factorization scale dependence of the cross 
section, and jiu is the moment-space anomalous dimension of the MS density We have 

7i/j = — Ajln Ni + 7, with 7, the parton anomalous dimensions. 

One can evaluate the anomalous dimensions and matching functions appearing in the re- 
summed cross section (1271) order-by-order in perturbation theory. The NNLL calculations 
mentioned at the beginning of the section require the anomalous dimensions A at three loops, 
all other anomalous dimensions at two loops, and the hard and soft functions at NLO. Obvi- 
ously, every time a term is added to the anomalous dimensions a whole tower of logarithms in 
the fixed-order expansion is resummed into the exponent, which accounts for the nomenclature 
behind this re-organization of the perturbative series. A feature of the resummed cross section 
in Mellin space is that it contains factors of a s evaluated below the QCD scale Aqcd, and 
is thus subject to Landau-pole ambiguities related to how to deal with this singularity. We 
will not discuss this issue in detail, but one way of circumventing this problem is to instead 
construct fixed-order expansions of fl2"7j) . a topic we return to in Section I3T31 
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3.2 Resummation with SCET 



We will explain the SCET approach to resummation using PIM kinematics as an illustrative 
example [13]. The same techniques apply to 1PI kinematics after only small modifications 
|69j . The starting point is the observation that partonic cross sections near threshold factorize 
into a product of a hard function, related to virtual corrections, and a soft function, related 
to real emission in the soft limit [9]. We write this factorization, valid up to corrections of 
0(1 - z), afl 



d&jj (z, Ma, m t , cos 6, fi F ) 
dMftd cos 9 



Tr 



HijiMtt, m t , cos 9, /i F ) Sij(VI(l 



z), m t , cos 9, fi F ) . (29) 



The boldface indicates that the hard functions Hij and soft functions Sij are matrices in color 
space. In SCET, the separation of the partonic cross sections into a product of hard and soft 
functions is achieved by a two-step matching procedure from QCD to the effective theory. Hard 
fluctuations are integrated out in a first step and the hard matching functions are identified 
as Wilson coefficients of an operator built of soft and collinear fields. Soft fluctuations are 
integrated out in a second step and the soft matching functions are the Wilson coefficient of a 
collinear operator, whose matrix element defines the PDFs. Resummation of large logarithms 
is then performed by deriving and solving the RG equations for the two functions. The main 
technical complication is that the RG equation for the soft function is non-local. This is 
dealt with using the Laplace-transform technique introduced in [70], which observed that the 
Laplace transformed functions 



s(L, M, m t , cos 9, a) = —= [ du exp f -—- J S(oj, M, m t , cos 9, u) 



(30) 



can be shown to obey local evolution equations whose solution can be inverted back to mo- 
mentum space analytically. 

The final expression for the resummed partonic cross section in PIM kinematics is 



da(z, M t t, m t , cos 9, fx F ) 



dMtjd cos 9 



exp 



Ml 



x Tr 



U(M, m t , cos 9, fi h , (i s ) H(M, m t , cos 9, fi h ) Lf\M, m t , cos 9, fi h , /i s 



x s (d v , M, mt, cos 9, ji s 



T(2 V ) (1-z) 



2E 



PIM/ 



2'/ 



ri=2a T ([i s ,[i F ) 



(31) 



In addition to the hard and soft matching functions, the formula contains factors related to 
the RG evolution from the matching scales fih and fi s to the factorization scale fj,p. These RG 
factors are given in terms of integrals over anomalous dimensions; their exact definitions can 



1 The parametric scaling M t i ~ m t is assumed , which is valid as long as the top quarks are not too highly 
boosted. 
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be found in jT3j. For values fi s < \ip the parameter r\ < 0, and one must use a subtraction 
at z = 1 and analytic continuation to express integrals over z in terms of plus distributions. 
Formula (13 ip can be evaluated order-by-order in RG-improved perturbation theory, using the 
standard counting ln/x/j//i s ~ ln(l — z) ~ l/« s . The current state of the art is NNLL |13j . 
which roughly speaking requires the soft anomalous dimension at two loops (as obtained in 
f27\ 12"%] ) and the hard and soft matching functions to NLO. 

Note that the logarithmic corrections at the soft scale are generated by derivatives with 
respect to r\ acting on the factors of (2Ef lM (z) / ii s ) 2rt ■ In the SCET calculation, it is natural 
to use 2Ef 1M (z) = M t t(l — z)/y/~z, which is the energy of the soft radiation in the partonic 
center-of-mass frame. We will refer to such a choice as the PIMscet scheme. An equally 
appropriate method, taken in all earlier literature on soft-gluon resummation, is to use the 
expansion of the energy in the z — > 1 limit, and we will refer to such a choice as the PIM 
scheme-the situation is summarized in Table EJ The difference between the two schemes 
involves power-suppressed terms of the form ln n z/(l — z) when the formula is re-expanded in 
fixed order. Such terms appear indeed naturally in the fixed-order calculations, as observed 
in the case of SCET applications to Drell-Yan [7T] and Higgs production [721 [75] . At NLO it 
can be shown numerically that keeping such terms in the factorized cross section (|29l) reduces 
the size of the remaining power-suppressed corrections which vanish in the z — > 1 limit. We 
will come back to this point in the discussion of approximate NNLO formulas in Section 13.31 

The methods of factorization and resummation described above for the case of PIM kine- 
matics generalizes to that of 1PI kinematics in a straightforward way. In fact, the structure 
of the resummed formula is exactly as in (13~T|) . The main difference is that the expression for 
the soft function changes, since it is derived by integrating over a different region of the phase 
space. The hard function, on the other hand, is the same, as is the form of the evolution 
equations. An analysis to NNLL, which required the calculation of the NLO corrections to the 
soft function, was performed in [T5]. Similarly to the case of PIM kinematics, logarithmic cor- 
rections in 1PI kinematics are generated by derivatives acting on a factor of (2E\ m (s±) / H> s ) 211 
and one has the choice between identifying this factor with the exact energy of soft-gluon 
radiation in the t rest frame or its expansion in the S4 — > limit. The two schemes defined 
in this way are summarized in Table [2j The difference between the IPIscet and 1PI schemes 
involves terms of the form ln"(l + s^/m 2 )/s^ when expanded in fixed order. Even more so 
than for the case of PIM kinematics, these subleading terms make an important numerical 
difference. 

The SCET approach to soft-gluon resummation can also be applied to the total cross 
section in the /3 — > limit. In this limit one can perform soft gluon resummation through 
the same techniques as in PIM and 1PI kinematics. However, such an approach is insufficient 
because one must also take into account Coulomb singularities of the form ln m /3//3 n , and a 
joint resummation of the two effects is more complicated than soft gluon resummation alone. 
This issue was solved in [Ml EI], which used techniques in SCET and non-relativistic QCD to 
perform a combined resummation of Coulomb and soft-gluon effects to NNLL order. 

From a technical standpoint, resummation within the framework of SCET is very similar 
to the Mellin-space resummation described in the previous section. In fact, the resummed 
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PIMsCET 


2Ef M {z) 


= M a {l - z)lyfz 


IPIsCET 


2Ef l 


>4) 


= s 4 /^mf + s 4 


PIM 


2Ef 1M {z) 


= Mu{l - z) 


1PI 


2El P1 


>4) 


= s 4 /m t 



Table 2: The values of the parameter E s which define different calculational schemes used in 
this review. 

formulas are completely equivalent if consistently re-expanded to any given order in a s 1 Other 
than the fact that the SCET results apply directly in momentum space, the main difference 
between the two approaches is the way in which the soft matching scale [i s is chosen, and the 
reasoning underlying this choice. In the SCET approach, one argues that the appearance of 
a well-separated soft scale at the level of hadronic cross sections is a dynamical effect due to 
the sharp fall-off of the PDFs away from the partonic threshold region [71]. To determine its 
numerical value, one studies the corrections arising from the soft function to the differential 
hadronic cross section as a function of fi s in fixed order, and finds the value at which these 
corrections are minimized. This numerical value is a function of the kinematic variables 
observed in the differential cross section (for instance in PIM kinematics of the invariant mass 
M t t), and is interpreted at the scale at which the soft function is free of large logarithmic 
corrections. Such a choice eliminates the Landau-pole ambiguity, but also implies that the 
SCET resummation exponentiates a subset of the higher-order plus distributions along with 
logarithms of the numerical ratio \i s j\ip\ a more complete discussion can be found in the 
Appendix of [T5] . 

3.3 Approximate NNLO formulas 

In the previous two sections we focused on all-orders resummation formulas. The basic idea of 
such formulas is to use properties of real radiation in the soft limit to calculate an infinite set 
of logarithmic corrections to the partonic cross section in terms of a smaller group of objects 
such as anomalous dimensions and matching functions. An alternative is to use the same 
formalism to determine only the logarithmic corrections to a certain accuracy in the fixed- 
order expansion. This is a useful approximation if it can be argued that such logarithmic terms 
capture the dominant corrections at a given order, and also that the higher-order logarithmic 
corrections are not so large as to spoil the convergence of the fixed-order expansion. A benefit 
of using such an expansion is that the Landau-pole singularities in the Mellin-space formalism 
are absent, and there is no need to introduce a dynamically generated numerical soft scale in 
the SCET approach. In this section we explain the structure of approximate NNLO formulas 
derived from NNLL resummation within the different types of soft limits. We point out that 
such formulas are subject to a number of ambiguities in the treatment of subleading terms in 
the soft limit, and discuss how these are dealt with in the literature. 

We begin with a discussion of 1PI and PIM kinematics. The approximate NNLO formulas 
in these cases are applicable to the differential distributions shown in Table [TJ but to compare 

2 This was shown in detail in [70l[74], and later in many other cases, including heavy-quark production in 
the production threshold limit |5T]. The same techniques apply to PIM and 1PI kinematics but an analysis 
has not yet appeared in the literature. 
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with the production threshold limit we will work at the level of the total partonic cross section. 
In PIM kinematics, one can write the general form of the NNLO corrections to the partonic 
scaling functions (fl3"]l as (suppressing the labels for the qq and gg channels) 



d cos Odz 



PIM 



n=0 



ln n (l - z) 
l-z 



C plM 5(l-z) + R^ M (z) 



PIM, 



(32) 



while in IPI kinematics, one can write 



/ (2) (/3,a0 = / diids 



n=0 



ln n (s 4 /m^ 
s 4 



+ C lFl 6(s A ) + R lFl (s A ) 



(33) 



The Di, C, and R coefficients are functions of the variables s, ti, Ui, m t , fi, so one evaluates 
the formulas above after an appropriate change of variables. The definition is such that the R 
coefficients are regular in the limit z — > 1 or s 4 — > 0. Soft gluon resummation at NNLL accuracy 
determines all of the Di coefficients in the limit z — > 1 or s 4 — > 0, and also the /i-dependent 
piece of the delta- function coefficient C [T4"l [T5l 175] . To determine the /i-independent piece of 
the 5-function coefficient would require the calculation of the soft and hard matching functions 
at NNLO order. The regular coefficients R are not determined by the soft gluon resummation 
in the form discussed here and can only be obtained through a complete NNLO calculation. 

An alternative to the IPI and PIM results is to approximate the NNLO corrections through 
the singular terms in the production threshold limit — > 0. The complete answer was first 
obtained in [76]. To illustrate its structure, we decompose the NNLO correction to the scaling 
functions (TT31) as 



4? = /T + /Ti» 



/4 



f, 



(2 ' 2) In 2 
ij 



/4 

2 

mr 



(34) 



Up to pieces which are regular in the /3 — > limit, the scale- independent corrections read [TB] 



(2,0) 



%(3 



qq 



3 60774 1 

— — — + - f - 140.368 In 2 /3 + 32.106 ln/3 + 3.95105 



(16vr 2 ) 2 9 

+ 910.222 In 4 (3 - 1315.53 In 3 (3 + 592.292 In 2 (3 + 528.557 In f3 



f(2,0) 

J t 



1 



77T/3 



1 + 

r 



(16tt 2 ) 2 192 L f3 2 
+ 4608 In 4 (3 - 1894.91 In 3 (3 - 912.349 In 2 (3 + 2456.74 In (3 



+ 



(35) 



Earlier results, which give very similar numerical results for the total cross section although the 
analytic structure is incorrect, were obtained in [51]. In that paper the coefficients multiplying 
the //-dependent logarithms in fl34")) . exact in (3, were also presented. 

In Section H] we will compare results for the total top-pair production cross section based 
on three approximate NNLO calculations available in the literature: the IPI threshold as 
implemented in [H], the IPI and PIM thresholds as implemented in the computer program 
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from [15], and the production threshold as implemented in the HATHOR program [77]. Since 
[H] and [15] both use 1PI kinematics and therefore differ only through implementations, it 
is important to clarify the sources of difference between the 1PI results from [H] and those 
from the SCET-based approach in [15J. Here there are several issues to consider. First, the 
numerical results depend on how the relation s + U\ + 1\ = is used in the plus-distribution 
coefficients (1331) before numerical integration, and the two papers in general use different re- 
writings which are equivalent in the limit —> 0. Second, even if this relation were used in 
the same way, there would still be differences in the analytic structure of the approximate 
NNLO formulas. These definitely affect the /^-dependent part of the 5-function terms, but 
possibly also the Do an d /i-independent part of the (^-function terms. The expressions for these 
coefficients used in [H] can be obtained from [TUl [TBI EH]; while those from [T5] were specified 
in that paper and given in the form of a computer program. Finally, the authors of [15] prefer 
the IPIscet (and PIMscet) implementation of the approximate NNLO formulas (see Table |2] 
and the discussion in Section l3~2]) . where a certain set of subleading terms appearing naturally 
within the SCET formalism is kept in the regular coefficients R. On the other hand, [2] uses 
damping factors 2m t / V§ in the NNLO soft-gluon corrections to the total cross section (for the 
Pt distributions a damping factor 2m±/\/§ is used) in order to reduce the contribution from 
kinematical regions further away from threshold, and thus improve the 1PI approximation. For 
the total cross section, the damping factors lead to differences between [H] and [15] that are 
actually not subleading in the S4 — > limit of the NNLL calculation. At NLO, the agreement 
between the exact result for the total cross section and its threshold approximation is improved 
in both the IPIscet scheme and the 1PI scheme with damping factors compared to the pure 
1PI expansion. However, at NNLO these different schemes produce rather different results, 
and the numerical effect of using damping factors versus including subleading SCET terms 
is greater than that from the other differences in implementation mentioned earlier in this 
paragraph. More details can be found in the numerical studies in [2] and [To], and also in 
the comparison in the next section. 

4 Total top-pair production cross section 

The most basic quantity in top-quark pair production is the total inclusive cross section. In this 
section we collect numerical results for this quantity based on different levels of perturbative 
accuracy and discuss the associated uncertainties. We focus on a comparison of NLO and 
approximate NNLO predictions derived from NNLL resummation in three types of soft limits: 
production threshold, PIM, and 1PI. 

To provide a QCD prediction for the total cross section at a particular collider, one must 
decide on two things. First, the value of input parameters such as the top-quark mass and a 
PDF set, and second, a procedure for estimating the theoretical uncertainties associated with 
the perturbative corrections beyond the accuracy of the calculation. For the NLO calcula- 
tion, one typically uses scale variations to estimate the uncertainties related to uncalculated 
corrections at NNLO and beyond. For the approximate NNLO calculations the uncalculated 
terms include the pieces of the NNLO correction which are subleading in the limit in which the 
formula is derived. The situation is thus slightly different from a full fixed-order calculation 
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and there is no conventional way for estimating theoretical uncertainties: some authors use 
the standard method of scale variations, some use a more complicated procedure. We describe 
the different approaches in conjunction with the results below. 

The results for the cross section in the pole scheme are summarized in Table |3j In addition 
to the NLO results, we use the following implementation of approximate NNLO predictions: 

• Production threshold results as obtained by the HATHOR program [77J with the default 
settings. 

• 1PI results as obtained in [T3] , 

• IPIscet and PIMscet results as obtained in [15J. These are combined into a final result 
for the cross section using the procedure and computer program presented in [69] ■ 

By default, we have set \i = = /ip = m t , with m t = 173 GeV. We display NLO results 
using MSTW2008 NLO PDFs, while for approximate NNLO results we use MSTW2008 NNLO 
PDFs [79]. Uncertainties in the HATHOR and 1PI results from [T3] are estimated by varying 
ji up and down by a factor of twcll, while uncertainties in [69] are estimated by independent 
variations of /ir and [if by factors of two, along with a scan over the values of the cross section 
in PIM and 1PI kinematics. 





Tevatron 


LHC (7TeV) 


NLO 
Aliev et. al. [77] 
Kidonakis [H] 
Ahrens et. al. [69] 


a 7/1+0.36+0.37 
u - ' ^-0.76-0.24 

7 1 q+0.31+0.36 
' • ±d -0.39-0.26 

7 no+0.00+0.36 
' - UO -0.24-0.24 

R CC+0.08+0.33 
U.UU_q 41 _ .24 


i6o±i;±i 

163™ 

156±|±« 



Table 3: Results for the total cross section in pb at NLO and within the various NNLO 
approximations. The first uncertainty is related to perturbative uncertainties, and the second 
is the PDF error using the MSTW2008 PDF sets [79] at 90% CL. 

An examination of the numbers in the table reveals the following features. First, the 
perturbative uncertainties in the NLO result are on the order of 20% at both the Tevatron 
and the LHC. This is a bit larger than the PDF uncertainty in both cases, although especially 
at the LHC one may obtain rather different results with other PDF sets, we refer the reader 
to [HO] for a recent discussion of this issue. Second, the perturbative uncertainties in the 
approximate NNLO results as obtained through the individual calculations are invariably 
smaller than at NLO-depending on the implementation, the uncertainties are reduced by a 
factor of roughly two to three, and are thus under the PDF uncertainties. At the LHC, the 
different NNLO approximations are in relatively good agreement, though the cross section of 
[69] is somewhat smaller than in [T^|[77]. At the Tevatron, on the other hand, the results from 
[HJ [77] are significantly larger than those from [69J. In fact, the range of values spanned by the 
three different approximate NNLO results at the Tevatron is about as large as that spanned 

3 Note that independent variations of \xf and /jr in fTJ] do not increase the uncertainty for LHC energies. 
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by the NLO calculation. Given the discrepancy, one is faced with the choice of estimating 
the theory uncertainties through the NLO calculation, with the spread of approximate NNLO 
values from the three different calculations, or by a particular NNLO approximation alone. 
The authors of [69| [TU [77] all give arguments in favor of their particular implementation of 
soft-gluon resummation, but it is beyond the scope of the review to properly summarize them. 
Instead, we describe very briefly the reasons behind the differences and explain the points 
which must be addressed to argue for one implementation over the others. 

The most important discriminator between the approaches is of course the soft limit in 
which the resummation is performed. As mentioned before, the production-threshold limit 
f3 — > is actually a special case of the 1PI and PIM thresholds. At the level of the approximate 
NNLO formulas, the production threshold results can be obtained from the 1PI and PIM 
results by re-expanding them in that limit, up to differences related to Coulomb terms (of 
the form ln m /3//3 n , see [13] for explicit results) which turn out to be very small numerically 
[TT| IT3] . particularly at the Tevatron. Therefore, in cases where the production threshold 
limit differs substantially from the 1PI and PIM cases, it can be concluded that subleading 
terms in /3 are not generically small, so arguing for the HATHOR results requires explaining 
why power corrections to the partonic cross sections do not inherit this generic property and 
furthermore why the subleading contributions included in the 1PI and PIM results are non- 
physical. Several differences between the 1PI results of [H] and the SCET-based results of 
[69] were described in Section 13.31 The numerical importance of various differences is not 
obvious from the table, since [HH] uses a scan over 1PI and PIM kinematics as an estimate of 
power corrections, but a more detailed study reveals that the dominant effects are the damping 
factors in [U] or the subleading terms in the IPIscet scheme in [69] , and less so the differences 
in the analytic structure described in Section 13.31 Favoring the results of either [H] or [69] 
thus requires arguments as to why the 1PI results with damping factors or the combination of 
the SCET implementations of the 1PI and PIM results gives a more reliable estimate of the 
cross section and its uncertainties. 

Recent measurements of the production cross section have been made at the Tevatron 
CDF [SI [821 [831 EH ES] and DO [23 E2 EE1 E21 EDI EI] experiments, and the LHC ATLAS 
[92~j 193"! IM] and CMS [951 ESI EZ] experiments. An especially useful method for comparing 
with the theory predictions is to show both as a function of the top-quark mass. Then the 
regions of overlap give a value of the top-quark mass as determined from the production cross 
section. While the errors in the top-quark mass extracted in this way are larger than those 
from kinematic distributions in top-quark decay, an advantage is that the theory calculations 
are carried out in a well-defined renormalization-scheme for the top-quark mass. The results 
of such analyses carried out through measurements of the lepton + jets channel by the DO 
collaboration at the Tevatron appear in [SHI Effl EB]- An analogous study has been performed 
by the ATLAS collaboration at the LHC in [99] . It is also interesting to perform the analysis 
using short-distance masses such as the MS [T2~l 169] or threshold masses [69]; results in the 
MS scheme extracted from Tevatron data can be found in [98] . In Figures [2] and [3] we com- 
pare recent experimental measurements assuming a pole mass of m t = 172.5 GeV with the 
NLO and approximate NNLO predictions, using the same method of estimating perturbative 
uncertainties as described above and used in Table [3J Although at the Tevatron the central 
values for an assumed top-quark mass of 172.5 GeV favor the higher values of the cross sec- 
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Figure 2: Experimental and theoretical values of the tt production cross section in pb at 
the Tevatron, assuming a top-quark mass of 172.5 GeV. The central values are indicated 
by the dots, the horizontal lines in the experimental results represent uncertainties from a 
combination of statistical and systematic errors, and the horizontal lines in the theory results 
are the perturbative and PDF uncertainties added in quadrature. 



tion predicted by [HI [77], the lower values of m t extracted through the NLO and approximate 
NNLO calculations from [BH] are still within the world average [50] . 

5 Top quark differential distributions 

We now turn to calculations of differential cross sections. To keep in line with the spirit of 
this review we cover only a few cases: the tt invariant mass distribution and the top-quark p? 
and rapidity distributions. Results for the pair invariant mass distribution using soft gluon 
resummation at NLO+NNLL or approximate NNLO exist from the calculations of [13] within 
PIM kinematics, while those for the p? and rapidity distributions from the calculations of 
PUDS US] wi thin 1PI kinematics. 

We first discuss the pair invariant-mass distribution. This is an important quantity for phe- 
nomenology. It is sensitive to searches for new physics through narrow s-channel resonances, 
which would show up as bumps in the distribution, while within the SM itself moments of 
the distribution give additional information about the top-quark mass [101] . Moreover, recent 
measurements of the forward-backward asymmetry by the CDF collaboration at high values 
of the pair-invariant mass are much larger than SM predictions, and an important constraint 
on new physics models which could produce such an excess is the level of agreement of the 
invariant-mass distribution with the SM calculations. At the Tevatron, this level of agree- 
ment is quite good, which we illustrate through the comparison of the CDF result with the 
NLO+NNLL calculation [13] in Figure [U Obviously, there are no narrow resonances in the 
experimental distribution at the Tevatron, and measurements at the LHC will extend this 
search to higher values of the invariant mass. It is worth emphasizing that at very large val- 
ues of the invariant mass the theory calculation is more complicated. This is mainly because 
the top quarks are highly boosted at very high invariant mass, and at some point one must 
apply an effective theory appropriate for the parametric scaling m t <C M t i. Also, for very 
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Figure 3: Experimental and theoretical values of the ti production cross section in pb at the 
LHC with y/s = 7 TeV, assuming a top-quark mass of 172.5 GeV. The central values are 
indicated by the dots, the horizontal lines in the experimental results represent uncertainties 
from a combination of statistical and systematic errors, and the horizontal lines in the theory 
results are the perturbative and PDF uncertainties added in quadrature. 



precise predictions, one must account for electroweak corrections, since these contain Sudakov 
logarithms that grow with the invariant mass |102[ I1U3] . 

Next, we discuss the top-quark transverse-momentum distribution. Results for the pt 
distributions at approximate NNLO at the Tevatron and the LHC with two collider energies 
are shown in Figure [5], based on the calculations of [H]. These results use MSTW2008 NNLO 
PDFs with fip = fiR = m t = 173 GeV. The DO collaboration in Ref. [104] has published 
measurements of the top quark p? distribution and compared with various theoretical results 
assuming m t = 170 GeV. We show a comparison of the DO measurements with the approximate 
NNLO results from [T3] in Figure El using scale variation around the central value /xr = /xp = 
m t = 170 GeV, and MSTW2008 NNLO PDFs. The agreement between theory [U] and 
experiment is very good. In Figure [7| we compare the DO measurements with the NLO and 
NLO+NNLL calculations from [T5]. MSTW2008 PDFs are used, and the theory errors are 
estimated with scale variations about the default values /ir = /ip = 2m t , with m t = 170 GeV. 
The NLO+NNLL predictions from [15] are also in good agreement with the measurements, 
although the optimal scale-setting procedure for the resummed calculation at higher values of 
p T where the distribution is small was left as an open point in that work. 

Finally, the NNLO approximate top-quark rapidity distribution at LHC and Tevatron 
energies from [H] is shown in Fig. El using MSTW2008 NNLO PDFs with fi F = /j, R = 
m t = 173 GeV. The rapidity distribution is important for the top quark forward-backward 
asymmetry which we discuss in the next section. 

Many more plots showing theoretical results for the pt and rapidity distributions can be 

found in m nans]. 
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Figure 4: Comparison of NLO+NNLL results [J3] for the invariant mass spectrum with CDF 
measurements [100] . 
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Figure 5: The top quark transverse momentum distribution at the LHC and the Tevatron 
from |14j . 
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Top quark p T at Tevatron S =1.96 TeV m=170 GeV 
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Figure 6: The top quark transverse momentum distribution at the Tevatron from [TJ] with 
m t = 170 GeV and \i = m t /2,m t ,2m t compared with DO measurements [104] . 
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Figure 7: Comparison of NLO+NNLL results [15] for the transverse-momentum distributions 
with DO measurements |104j . 
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Top quark rapidity distribution NNLO approx m=173 GeV 




Figure 8: The top quark rapidity distribution at the LHC and the Tevatron from [To] . 



6 Forward- backward asymmetry 

Closely related to differential distributions is the top-pair forward-backward (FB) asymmetry. 
Its basic definition is 

pi N(yi > 0) - N(yl < 0) = a(y t > 0) - a(y t < 0) a\ 
FB N(yi>0) + N(yi<0) a(y t > 0) + a(y t < 0) 4' 1 ] 

where N is the number of events and y\ is the top-quark rapidity in Lorentz frame i. The defi- 
nition in terms of the number of events is convenient for experimental measurements, whereas 
in theory calculations one obtains the forward-backward asymmetric (c?a) and symmetric (erg) 
cross sections by integrating rapidity distributions over the appropriate values. Since each of 
these are obtained as a series in a s , one must further decide whether to re-expand the ratio of 
asymmetric to symmetric cross sections in the computation of Ap^, a point we will come back 
to in a moment. The basic definition above can also be modified to include cuts on kinematic 
variables, for instance by imposing that M t i or \y t \ is greater or less than a certain value. 

At a pp collider such as the Tevatron, charge conjugation invariance in QCD implies that 
N(yl > 0) = N(y l t < 0), so the FB asymmetry is equivalent to the charge asymmetry. At 
a pp collider such as the LHC, rapidity distributions are symmetric and the FB asymmetry 
vanishes. On the other hand, the production rates for top and antitop quarks at a given 
rapidity are in general different, so partially integrated charge asymmetries at the LHC do not 
necessarily vanish. Moreover, one can define non-vanishing total asymmetries with respect to 
variables such as \rjt\ — \t]t\, with t] the pseudorapidity, or \y t \ — \yt\- Various definitions and 
results for (partially integrated) charge asymmetries at the LHC can be found elsewhere in 
this volume |lU5j and will not be covered here. From now on, we will focus solely on the FB 
asymmetry at the Tevatron. 
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Experimental measurements of the FB asymmetry at the Tevatron have been made in 
the ti and pp (laboratory) framejf] The measurements of the total asymmetry obtained by 
the CDF collaboration using 5.3 fb" 1 of data are [106j A^ B = (15.0 ± 5.5)% (pp frame) and 
Ap B = (15.8 ± 7.5)% (ti frame). The quoted uncertainties are derived from a combination 
of statistical and systematic errors. The measurements are roughly 2a (la) higher than 
theoretical results in the laboratory (ti) frame reviewed below. Even more interesting is the 
measurement as a function of the top-pair invariant mass. After grouping the events in two 
bins corresponding to M t i < 450 GeV and M t i > 450 GeV, the asymmetry in the latter bin 
was measured to be Ap B (M t i > 450 GeV) = (47.5 ± 11.4)% which is about 3a higher than the 
theory predictions in the SM. Attempts to explain the discrepancies of the CDF results with 
the SM predictions are reviewed in |105j . Here, we present the most up-to-date results for 
the FB asymmetry within the SM itself, explaining what is known in fixed-order perturbation 
theory and within different frameworks for soft-gluon resummation. 

In the SM the FB asymmetry is due mainly to QCD effects. In fixed-order perturbation 
theory the asymmetric cross section first appears at 0(a^), in other words at NLO compared 
to the symmetric cross section. The contributions can be traced to certain types of diagrams 
in the qq channel and qg channels, explicitly identified and calculated in |107[ 1108] . The gg 
channel does not contribute to the asymmetric cross section at any order in perturbation 
theory, because it involves a symmetric initial state. If the ratio o^/Vs i n (J36]) is consistently 
expanded in a s , as we will assume to be the case unless otherwise indicated, then the Afb ~ a s 
at the first non- vanishing order. We will label this leading contribution as NLO, with reference 
to the order at which the differential cross section is needed. The total FB asymmetry at NLO 
is small in the SM, about 5% in the pp frame and 7% in the ti frame-exact numbers will be 
given below. The qq channel gives the dominant contribution, while that from the qg channel 
is much smaller numerically. 

Given the potential deviations between the measurements and the SM, it is especially 
important to estimate the effects of higher-order QCD corrections on the FB asymmetry. 
This can be done using the methods of soft-gluon resummation. We show results for the total 
FB asymmetry in the ti and pp frames at NLO and with soft-gluon resummation included in 
Table H We use MSTW2008 PDFs, with m t = 173.1 GeV for the calculations at NLO+NNLL 
and m t = 173 GeV at approximate NNLO in the pp frame (the numerical difference from 
the mismatch in the top-quark mass is negligible). The central values refer to the choice 

— /J>f — f^t, and perturbative errors in fixed order are estimated with correlated variations 
by factors of two; in the resummed result also the matching scales /i^ and fi s are also varied 
as described in |109j . In all cases the theory error is much smaller than the experimental one 
quoted above, and the main qualitative finding is that resummation has only a rather mild 
effect on the FB asymmetry. This is true whether one uses NLO+NNLL, or the approximate 
NNLO calculations of [T3~| H~5| H~6]FI In fact, the same conclusions hold for the invariant-mass 
dependent CDF analysis, where the asymmetry is studied in a bin below and above a cut at 
M t i — 450 GeV. The studies in |109] showed that such a binning roughly divides the total 

4 Experiment actually measures the asymmetry with respect to the rapidity difference yt — yt, but this is 
equivalent to (jSHJ) in the ti frame so we do not distinguish this as a separate observable. 

5 In comparing results it is important to note that [16] does not expand the FB asymmetry in a s , evaluating 
instead the denominator of ([56]) numerically at NNLO with no further expansion. 
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[%] 


NLO [109] 


7 qo+0. 69+0.18 
1 ■<J^_o,59-0.19 


A Ol +0.45+0.13 
4: - o - L -0.39-0.13 


NLO+NNLL [109J 


7 0^+1.04+0.20 
' '^-0. 67-0. 27 


ZL SS+0.20+0.17 
^• oo -0.23-0.18 


NNLO approx . [16] 




e 9+0.0 

J - z -0.6 



Table 4: The FB asymmetry in the it and rest frames. The first error refers to perturba- 
tive uncertainties estimated through scale variations, the second (where applicable) to PDF 
uncertainties. 

asymmetric cross section in half, and that the effect of resummation in each bin is roughly 
the same as for the results shown in the table here-in other words, it is a mild effect. Similar 
conclusions can be drawn from the NLO+NLL study of [110j . 

The moderate size of the higher-order QCD corrections as estimated by soft gluon resum- 
mation makes it important to consider other sources of uncertainty in the SM calculation. 
First, there are of course sub leading terms in the soft limit to which soft-gluon resummation 
has no access. The magnitude of these terms is estimated through scale variation but will 
be known for sure only once the fixed-order calculation is done at the next order. Second, 
electroweak corrections (EW) can also be important for the FB asymmetry. These were con- 
sidered already in |lU8j . where they were estimated to increase the asymmetry roughly by a 
factor of 1.1. Refinements related to additional contributions were made in [lllj . where a 
slightly smaller contribution was found, but very recently the calculations of |112] quote an 
increase of about 1.2 due to EW corrections. The corrections for binned analyses in M t i and 
pair rapidity difference were found to be roughly the same in the bins studied by experiment. 
These corrections are comparable or larger than the QCD uncertainties and thus cannot be 
ignored-however, they do little to change the 3a discrepancy between the SM and the CDF 
measurement in the high invariant- mass bin. 

7 Single top production 

The observation of single top quark production at the Tevatron |113[ 1114] 1115] 1116] 1117] 1118] 
1119] and at the LHC |120] 1121] has increased the need for reliable theoretical calculations of 
the cross sections for the corresponding production processes. The single top cross section is 
less than half of that for ti production but the backgrounds are considerable thus making the 
observation of single tops quite challenging. 

Single top quarks can be produced through three distinct partonic processes. One of them 
is the t-channel process that proceeds via the exchange of a space-like W boson, a second is 
the s-channel process that proceeds via the exchange of a time-like W boson, and a third is 
associated tW production. The leading-order diagrams are shown in Fig. 

7.1 ^-channel production 

The t-channel partonic processes are qb — > q't and qb — > q't. At both the LHC and the 
Tevatron the t-channel is numerically dominant. Calculations of NLO corrections for t-channel 
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Figure 9: Leading-order t-channel (left), s-channel, and tW (right two) diagrams for single 
top quark production. 



LHC 7 TeV 


t 


t 


t-channei |122j 
s-channel [123J 
tW [124] 


41.7io;2 ± 0.8 
3.17±0.06ijj;ljj 
7.8 ± 0.2+°| 


22.5 ± 0.5±g;2 
1.42±0.0llg;g? 
7.8 ± 0.2+°j 



Table 5: Results for single-top and single-antitop approximate NNLO cross sections at the 
LHC with m t = 173 GeV in pb. The first uncertainty is from scale variation and the second 
is the PDF error using the MSTW2008 NNLO PDF sets at 90% CL [79]. 

production at the differential level have been known for some time |125] (see also further recent 
NLO studies in [T2EJ E27J EH CE2HJ HHP] ). 

Theoretical calculations for single top quark production beyond NLO that include higher- 
order corrections from soft-gluon resummation appeared in |131[ 1132] at NLL and in |122] at 
NNLL. The NNLL theoretical expressions in |122] were used to derive approximate NNLO 
cross sections for t-channel single top or single antitop production at the Tevatron and the 
LHC. 

In Figure [TO] the t-channel approximate NNLO cross section from NNLL resummation is 
plotted versus collider energy, together with recent measurements from the Tevatron |116[ 1118] 
and LHC experiments |120[ 1121] . For a top quark mass of 173 GeV the i-channel single top 
quark NNLO approximate cross section at the Tevatron is 

al op ch (m t = 173 GeV, V* = 1-96 TeV) = 1.04±g;gg ± 0.06 pb (37) 

where the first uncertainty is from scale variation between m t /2 and 2m t and the second is 
the PDF uncertainty, calculated using the MSTW2008 NNLO PDF sets [79] at 90% CL. We 
note that the results for single antitop production at the Tevatron are identical to those for 
single top. 

For t-channel production at the LHC we note that the single top cross section is different 
from that for single antitop production. For m t = 173 GeV we list the results for LHC 7 TeV 
energy in Table 

7.2 s-channel production 

The lowest-order processes in the s-channel are of the form qq' — > bt, which include the 
dominant process ud — > bt as well as processes involving the charm quark and Cabibbo- 
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t-channel single top + single antitop cross section 
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Figure 10: The single top plus single antitop t-channel cross section versus collider energy. 



suppressed contributions. The QCD corrections for s-channel production at next-to-leading 
order (NLO) are known at the differential level |125] and they substantially increase the cross 
section and stabilize the dependence on the factorization scale. 

Approximate NNLO calculations from NLL soft-gluon resummation appeared in |13![ 1132] 
and from NNLL resummation in |123j . The soft-gluon corrections dominate the cross section 
and the NLO expansion of the resummed cross section approximates very well the complete 
NLO result for both Tevatron and LHC energies. Results based on SCET have appeared in 
|133j . In addition to the difference in formalism, Refs. |123] and |133] differ in their definitions 
of the S4 variable in the resummation; we refer the reader to the papers for a more detailed 
discussion. Here we just note that contrary to the study in |123] the NLO threshold expansion 
in the approach of |133] was found to poorly approximate the full corrections at LHC energies. 

For a top quark mass of 173 GeV, the s-channel single top cross section and its associated 
uncertainties at the Tevatron are |123] 



where the first uncertainty is from scale variation and the second is the PDF uncertainty at 
90% C.L. The results for single antitop production at the Tevatron are identical. 

The results for single top and single-antitop s-channel production at the LHC at 7 TeV 
with nit — 173 GeV are shown in Table [5j 

7.3 tW production 

The associated production of a top quark with a W boson, bg — > tW~, is sensitive to new 
physics and allows a direct measurement of the V t b CKM matrix element. This process is 




(38) 
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practically negligible at the Tevatron but it has the second highest cross section among single 
top processes at the LHC A similar process in physics beyond the Standard Model is associated 
production of a top quark with a charged Higgs boson, bg — > tH~ . The next-to-leading order 
(NLO) corrections to bg — > tW~ were calculated in |134j . 

NNLO soft-gluon corrections from NLL resummation were calculated in |131[I132] and from 
NNLL resummation in [124] . The NLO expansion of the resummed cross section approximates 
well the complete NLO result for both Tevatron and LHC energies. 

The results for tW~ and iW + at the LHC at 7 TeV with m t = 173 GeV are shown in 
Table [5j The NNLO approximate corrections increase the NLO cross section by ~ 8%. We 
note that the cross section for bg — > iW + is identical at both Tevatron and LHC to that for 
bg tW~. 



8 Conclusions 

We have reviewed QCD calculations for inclusive top-quark production at Tevatron and LHC 
energies. The total and differential cross sections have been known at NLO in QCD for many 
years, and more recently different implementations of soft-gluon resummation at NNLL order 
have also become available and applied to many observables. Compared to NLO calculations, 
phenomenological predictions based on resummation invariably show smaller dependence on 
factorization and renormalization scales, and produce moderate changes with respect to central 
values for a default scale choice. However, especially for the total pair-production cross section 
at the Tevatron, different forms of resummation can yield results which show only nominal 
agreement with one another within the quoted uncertainties (see Table [3]). We explained in 
detail the reasons for such differences at the level of approximate NNLO results, but at present 
the best way of implementing resummation for the total cross section is still subject to debate, 
as is the most reliable way of estimating the uncertainties. Ultimately, only a complete NNLO 
calculation will resolve this point. 
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